Guided Modes of Elliptical Metamaterial Waveguides 
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The propagation of guided electromagnetic waves in open elliptical metamaterial waveguide struc- 
tures is investigated. The waveguide contains a negative-index media core, where the permittivity, 
e and permeability fi are negative over a given bandwidth. The allowed mode spectrum for these 
structures is numerically calculated by solving a dispersion relation that is expressed in terms of 
Mathieu functions. By probing certain regions of parameter space, we find the possibility exists to 
have extremely localized waves that transmit along the surface of the waveguide. 

PACS numbers: 



Waveguides are structures that are typically designed 
to transmit energy along a specified trajectory with min- 
imal attenuation and signal distortion. When transmit- 
ting surface waves, this implies confining the traveling 
wave within or adjacent to the waveguide walls [llj. Vari- 
ous avenues can be pursued when attempting to improve 
on a waveguides' capabilities, including modifying the 
constitutive effective material parameters of the guide. 
The recent surge of interest in negative index of refrac- 
tion materials [3], or negative index media (NIM), has 
prompted a reanalysis of many conventional results for 
waveguide devices, and a subsequent search for exotic 
transmission characteristics when incorporating this par- 
ticular composite media, or metamaterial, into a variety 
of waveguide configurations. 
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A crucial feature of NfM, is the frequency dispersion 
of the permittivity, e, and permeability, /i, with e and fi 
simultaneously rendered negative over a particular band- 
width. This results in wave propagation in which the 
phase velocity and energy flow of an electromagnetic 
wave can be antiparallel, and also the possibility of a 
negative index of refraction, proposed long ago Q. The 
study of various NIM based open waveguide structures 
and resonators during the past few years has demon- 
strated a number of interesting effects: it was shown Q 
that a quasi 1-D bilayer resonator containing a NIM layer 
can be substantially smaller than the usual cavity size 
due to phase cancellation. For the case of a thin planar 
NIM waveguide, the TM mode can propagate for arbi- 
trary widths and posses a single mode for slow waves [1] . 
Guided TE modes were also shown to have electric field 
profiles containing nodes [5] , and exhibit a sign- varying 
Poynting vector [5|, |6| . Similar results were reported in 
circular NIM fibers [7 1 . 

Distortion of the circular dielectric waveguide into an 
elliptical guide, while maintaining the cross sectional 
area, has been shown to reduce attenuation of the domi- 
nant mode and modal degeneracy, allowing for prac- 
tical guiding of traveling electromagnetic waves. Atten- 
uation effects and power fiow expressions were achieved 
for wave propagation in a surface wave transmission line 
with an elliptical cross section [l^. It was found that 
some modes in the guide have lower attenuation than 
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FIG. 1: Cross section of the open elliptical waveguide struc- 
ture. Within the metamaterial core, ei and ni axe in general 
both negative, and frequency dispersive. The outer region is 
air; £2 = 1 and fj,2 — 1. 



the corresponding modes in a circular guide. The slow 
and fast hybrid mode spectrum in a metallic elliptical 
waveguide with a confocal dielectric lining was also cal- 
culated [l4j], demonstrating the potential for a surface 
wave transmission device. Moreover, the distribution of 
electromagnetic fields in a double layer elliptical waveg- 
uide was calculated to first order and revealed in some 
cases, dispersion solutions where the power and phase ve- 
locity directions are antiparallel[8|. Spurred by the recent 
theoretical and experimental advances in NIM waveg- 
uides and resonators, we examine in this paper, open 
NIM waveguides with elliptical cross sections, affording 
a greater flexibility in the parameter space determined by 
the geometry and material constraints. For a given ec- 
centricity, e, we calculate the permitted propagation con- 
stants over a range of frequencies. The allowed modes are 
separated into fast and slow propagating regions of the 
dispersion diagram, where wave localization is demon- 
strated in the form of electric field distributions. 

When searching for exact solutions to Maxwell's equa- 
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tions, it is convenient to work in a coordinate system in 
which the boundary of the structure coincides with one 
of the coordinates being held constant. It is thus ap- 
propriate for the geometry under consideration to work 
in an orthogonal elliptical coordinate system described 
by the coordinates ^ and rj, depicted in Fig. [TJ Ellip- 
tical coordinates are related to their rectangular coun- 
terparts via X = p cosh ^ cos r/ and y = psinh^sin?], for 
< ^ < CO and < 77 < 27r, where p = {a? — Ip-)^/'^ is 
the semifocal length of the ellipse expressed in terms of 
the semimajor and semiminor axes a and h respectively. 
The waveguide boundary is located at ^ = Coj and hence 
a = pcosh^o and h — psinhfo- It then follows that the 
eccentricity, e, is written e = [1 — (&/a)^]-^/^ = l/ cosh^o, 
such that < e < 1, with 6 = corresponding to a circu- 
lar cross section. In order to obtain the electric (E) and 
magnetic (H) fields, we must solve the vector Helmholtz 
wave equation obtained from Maxwell's equations. This 
yields three coupled second order differential equations 
for either E or H. The waveguide along the axial {z) 
direction is translationally invariant and thus the form 
of the equation governing the longitudinal Ez (or Hz) is 
identical to the scalar Helmholtz equation in elliptic co- 
ordinates, permitting a semi- analytic and more tractable 
solution. The fields are assumed to vary harmonically in 
time and propagation occurs in the positive z direction. 
The wave equation then reduces to the following form: 
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Ez^{tv)+k!Ez^i^,V)^0, (1) 



where the index i denotes the region, 
p-\/cosh(2^) — cos(2?7), fcf — [ii{uj)ei{u:)u? j c? — k\, and 
for the unattenuated modes of interest here, fc^ is the real 
propagation constant along the z-direction. We have also 
suppressed the usual Qy^\>(ikzZ — iojt) factor. A similar 
equation exists for Hz ■ Due to the cylindrical symmetry, 
the transverse field components (i?,,i, E^i) and {H^ii, H^i) 
can be determined [10] from Ezi and Hzf. 
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We proceed by inserting Ezi{£,,f]) — Ui{(,)Vi{r]) into 
Eq. ID), splitting it into two ordinary differential equa- 
tions with separation constant A: 



d^V^{r]) 



+ (A - 2qi cos2r|)V^{r|) = 0, (3a) 



-(A-2<7,cosh20C/.(0 =0, (3b) 



where qi = (fcip/2)^, and 92 = —(^2^/2)^. 

The angular Mathieu equation [Eq. (f3a|) ] describes the 
angular variation of the field around the ellipse. Two 
linearly independent periodic solutions exist, akin to 
the trigonometric sin and cos functions: the even and 
odd angular Mathieu functions, denoted cen{r];qi) and 
se„(?7; qi), respectively. For the elliptical waveguide prob- 
lem, the angular Mathieu functions must be periodic with 
period 2tt and of integer order, otherwise the solution set 
is nonperiodic and can be unstable [l^. We can thus 
expand the angular Mathieu functions in a Fourier se- 
ries, where the expansion coefficients can be found recur- 
sively by substituting the expansions back into Eq. (|3ap . 
The parameter A found in Eqs. ([^can be solved using 
a method of continued fractions [la]. The orthogonal- 
ity and normalization relations of the angular Mathieu 
functions are given by Eq. (|XT|) . The angular Mathieu 
functions are also periodic in 77 only for special charac- 
teristic values of A, denoted here {qi) for the even 
solutions and b„{qi) for the odd ones. 

The radial Mathieu equation (|3bp. admits two gen- 
eral types of solutions: the radial Mathieu functions of 
the first and second kind. The solutions of the first 
kind are divided into even and odd functions denoted 
as, Cen(£,',qi) and Sen{£,\qi) respectively. The calcula- 
tion of the characteristics numbers and the Fourier coef- 
ficients are exactly the same as for the angular functions. 
In identifying a convenient series expansion to compute 
the radial Mathieu functions, a number of techniques are 
applicable |ll|. We found, in agreement with past works, 
that the Bessel J„ product series are the most stable ba- 
sis of functions to use [12]. Not surprisingly, the Ce„ 
and Sen functions coalesce into J„, as the elliptical cross 
section degenerates to a circular one. The Mathieu func- 
tions of the second kind for qi> are denoted Feyn{£,', qi) 
(even), and Geyn{£,',qi) (odd), and are calculated simi- 
larly, except the series expansion involve products of both 
the Bessel Yn and J„ functions. The functions Feyn and 
Geyn are analogous to the Bessel functions of the second 
kind, Yn, in circular coordinates. When qi < 0, Feyn 
and Geyn transform into Fekn{^;qi) and Gekn{£,;qi) 
(analogous to the Bessel K functions), and they are 
related [13], e.g., for even order, 2Fek2n{^', ^qi) = 
(-1)" [-Fey2n{^ + iTr/2; q,) + iCe2n{^ + iTr/2; q,)]. 

For a general cylindrical open guide with circular cross 
section, pure TE or TM modes exist only for symmetrical 
electromagnetic fields, i.e., independent of the azimuthal 
angle. For an open elliptical waveguide however, the re- 
duction in symmetry forces the electromagnetic modes 
in a given region to be hybrid in that both the longi- 
tudinal fields, Ezi and Hzi exist simultaneously (i de- 
notes the region). Based on the discussion above, it is 
clear that field solutions in an elliptical domain are split 
into even and odd components: Ezi {^zf-^zi} and 
Hzi {Hzi, Hzi}- The procedure adopted here consists 
of writing H^^ {H°^) waves as products involving even 
(odd) Mathieu functions, and E°i {E°^) waves in terms of 
odd (even) Mathieu functions 9] . The complete solution 
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is then expanded as products of angular and radial Math- 
ieu functions of the requisite parity, and that obey the 
radiative condition. The boundary conditions constrain 
the types of Mathieu functions used in the expansions 
[l^ . as we require "stationary waves" in the transverse 
r] and ^ direction. Within each region, we thus seek so- 
lutions of the form, 



£;j,(C, r],z)=Y^ a,niArm{e, qi)se.m{m q^)e''''", (4a) 

m— 1 

oo 

Ht^it br,iiB,mie, g»)ce™(7?; g.)e*'''^ (4b) 
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where ami and b„ii are constants, and the coefficients 
Am,i{£,', Qi) , and Bmi(S,](li) contain the "radial" depen- 
dence to the fields, given below. Without loss of gen- 
erality, the methodology used here will focus on waves of 
even parity, as the procedure for odd waves is similar. 

Within the waveguide region (see Fig. [1]), the solu- 
tions are the radial Mathieu functions of the first kind, 
^mi(C;9i) = Sejn{^\qi), and B„ii(^;gi) = Cem{^\qi). 
For the guided modes of interest here, the fields must 
decay from the surface at ^ — ^o, and therefore 
within the surrounding medium we have, ^m2(C;92) = 
Gekmi^;-q2) and B™2(C; 92) = Fe/cm(^; -92)- With this 
requirement on the fields, and for a given set of mate- 
rial and geometrical parameters, a restricted number of 
waveguide modes exist. To determine the allowed modes, 
the tangential E and H fields are matched at the bound- 
ary ^ = ^0 separating the two media. The 77-dependence 
is integrated out by making use of the orthogonality prop- 
erties of the angular Mathieu functions (see Appendix) . 
We then cast the boundary matched equations into a lin- 
ear equation system containing in principle, an infinite 
hierarchy of Mathieu functions. The higher order Math- 
ieu functions arise from the lack of one-to-one correspon- 
dence between the angular Mathieu functions in regions 
with differing material parameters: the arguments of the 
angular Mathieu functions depends on the parameter qi, 
which in turn depends on e and fi of the relative media. 
This is in contrast to a circular waveguide, where the an- 
gular dependence is a function of only integer multiples 
of the azimuthal coordinate, m(j). 

In order to find the nontrivial solution, the problem of 
finding the allowed modes thus amounts to finding where 
the associated determinant of the linear equation system 
vanishes over a range of frequencies, propagation con- 
stants and cUipticities. The methodology we shall discuss 
is valid for hybrid HEn or EHu (the first letter repre- 
sents the dominant field) modes [9|, with small e, and 
frequencies corresponding to small |ei — £2! (and small 
|a*i ~ M2|)- Under these conditions, the expansions in 
Eqs. (m can be limited to the first few terms. Taking for 
example, the two lowest order terms in the outer region. 



yields the following even mode dispersion relation, 
I ei Se[{qi) S2 Gek[{-q2)\ I Hi Ce[{qi) Fek[{-q2) 
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where kg = lj/c, Ci = -/3ii7ii - Tis/^ai, and C2 = 
aii7ii+a3i73i- The quantities a^n, t™„, Vmn, and 
7m„, outlined in the Appendix, arise from the integrals 
of products of overlapping angular Mathieu functions. 
Writing the dispersion relation explicitly in this way re- 
duces computational time considerably by avoiding un- 
stable numerical determinants of large matrices and the 
associated multiple function calls to higher order Math- 
ieu functions. The odd mode spectrum is easily obtained 
via the interchange <-!■ — /i^. As the elliptical cross 
section degenerates into a circular one, Ci 1, C2 
— 1, and the Mathieu functions appropriately transform 
into their corresponding cylindrical Bessel functions: 
{Ce[{qi)/Cei{qi),Se[{qi)/Seiiqi)} ^ uJi{u)/Mu), 
and {Fek[{-q2)/Feki{-q2),Gek[i-q2)/Geki{-q2)} 
vK[{v)/ Ki{v). In this limit, the dispersion relation ([5]) 
reduces to the familiar characteristic equation for a cir- 
cular waveguide: 
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where — kfa^ and v"^ — — fc|a^. 

The surface wave dispersion relation ([5]) is a function 
of the parameters, k^, uj, and eccentricity, e; only partic- 
ular combinations of these quantities that satisfy Eq. ([5]) 
are allowed mode solutions. The hybrid waves that are 
explored here can possess even and odd components, 
and are denoted appropriately in subsequent figures. We 
present the propagation constant in terms of the conve- 
nient dimensionless ratio, k^/ko, and the frequency units 
are all in GHz. The waveguide cross section is assumed 
to not deviate greatly from that of a circular guide, re- 
fiected in moderate values of e. The permittivity and 
permeability in the NIM regions, Si and /ii, respectively, 
have the frequency dispersive form, ei = 1 — (ujp/uj) and 
/ii = 1 — (o^m/cj)^ , where Up and uim are the effective elec- 
trical and magnetic plasma frequencies respectively. 
Thus for oj < ujp/^/2, we have £1 < —1. The study here 
is concerned with frequency regions of parameter space 
where fj,i and £1 are simultaneously negative. 

Determining the allowed modes typically involves hold- 
ing the semimajor and semiminor axes of the waveguide 
(and hence e) fixed, and then scanning Eq. ^ over kz 
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FIG. 2: Surface wave spectrum for the even and odd hy- 
brid modes of an elliptical waveguide filled with NIM and 
embedded in air (top panel). The bottom panel is the same 
configuration with positive The propagation constant, kz, is 
normalized by the free space wavenumber fco and is plotted as 
a function of the operating frequency for a cross section corre- 
sponding to an ellipticity of e = 0.44. We take the character- 
istic frequencies in the material dispersion to be ujp — 8(27r) 
GHz, and u>m ~ 7(27r) GHz, which corresponds to ^i, ei < — 1 
within the frequency range, / < Lu„, /{2nV2). The mode 
spectra illustrates the slow and fast modes, separated by the 
gi = curve. If the NIM core is replaced with media corre- 
sponding to positive ei and /xi (bottom panel), the slow wave 
solution is absent and the curves vary in an opposite sense. 



and CO. Any sign change that occurs signifies a zero cross- 
ing that can be pinpointed through an iterative root find- 
ing process. Other variable combinations may be used, 
depending on the parameter study. In the top panel 
of Fig. [2] we show the mode spectra found by solving 
Eq. ([5|), and its odd counterpart, for a waveguide with 
eccentricity, e = 0.44. For clarity, only four neighbor- 
ing sets of dispersion curves within the given frequency 
window are shown. Two sets reside completely beneath 
the qi — curve, one set entirely outside of it, and an- 
other set that traverses both regions. For this geometry 
and range of frequencies, the paired even-odd mode solu- 
tions follow similar trends, separating at higher frequen- 
cies. The solutions to the Helmholtz equation (Eq.[T]) de- 
pend on the eccentricity, e, of the elliptical guide through 
the parameters qi and q2. We found as the ellipse flat- 
tens, more splitting occurs between the even and odd 
modes. In general, each branch of the even and odd 
modal curves coalesce at cutoff, where the propagation 
constant approaches the free space value {kz ko). The 
dispersion diagram portrays the allowed electromagnetic 
wave solutions that travel along the guide and it eluci- 
dates important information regarding the possible lo- 



calization characteristics of guided modes. In particu- 
lar, the first three sets of dispersion curves that satisfy 
^e2[ii < kz/ko < y/ei/j,i correspond to conventional 
surface waves. Within this parameter space region, the 
phase velocity of guided waves, Vp, exceeds the phase 
velocity of waves in a homogeneous bulk medium, i.e., 
Vp > c/^Sifii. These fast wave solutions to Maxwell's 
equations will decay in the air region, but not necessarily 
in the waveguide core. Indeed, as the structure increases 
in size, or as / decreases (increasing |ei|), the electric 
and magnetic fields inside the waveguide oscillate with 
shorter wavelengths. At cutoff, k2 ^ (thus q2 0), 
and the ensuing decay length increases outside the guid- 
ing surface. This causes a significant portion of the en- 
ergy flow to occur in the air region, where the flelds can 
then become less sensitive to the relevant geometrical pa- 
rameters, such as the ratio a/b. 

We see in Fig. 2, that for the case of a NIM core, 
there are also solutions that reside outside of the <7i = 
boundary. These slow wave solutions are noticeably ab- 
sent when Si and fii are strictly positive (see bottom 
panel). The presence of NIM sets up a non-oscillatory 
field profile that rapidly decays outside the guide, allow- 
ing for the possibility of guided modes that are more 
localized to the surface, akin to surface plasmon waves 
on metal surfaces studied long ago. Note that as k^ in- 
creases, 52 typically increases, in which case the radial 
Mathieu functions decay sharply, confining the field more 
to the waveguide. The slow wave solutions seen in the 
inset of Fig. [5] have dispersion characteristics that de- 
pend strongly on frequency. Near cutoff, these modes 
demonstrate that the direction of the group velocity, 
duj/dkz, evaluated at a particular fc^, can differ among 
the even and odd wave solutions. This is related to 
the (time-averaged) local energy-density flow along z, 
Szi = [l/'2.)^{E^.,H*i - Erj^H^i). An excited wave with 
the appropriate frequency, can reverse direction between 
layers, a hallmark of the peculiar wave propagation that 
can arise in NIM guiding structures. To address the re- 
lationship between net energy flow in the system and 
group velocity, we show in Fig. [3fa), the total power, P 
as a function of the calculated mode frequency, found by 
summing the power flow in the waveguide (Pi) and air 
(P2) regions. The power through a given cross sectional 
area. A, was calculated by integrating the z component 
of the Foynting vector over A, 
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1,2, 



(7) 



where h is the usual coordinate scale factor. The frequen- 
cies used in determining P are governed by the dispersion 
curve, shown in the inset of Fig. [3I^a). The arrows label 
points where dto/dkz is zero in the dispersion diagram, 
and correlate with zero net power flow in the system, 
i.e.. Pi = —P2- In general, we see from the figure that 
the direction of net power flow coincides with the sign of 
duj/dkz. A sign change in this slope causes the dispersion 
curves to bend back in the kz — lo plane, related to the 
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FIG. 3: (a) The total power, P = Pi + P2 (normalized by 
I Pi I + IP2I) versus the frequency for (odd) slow wave modes. 
The eccentricity, e, is set at e = 0.59. The propagation con- 
stant, fcz, increases in going from A — > B, in accordance with 
the mode dispersion diagram in the inset. In (b) energy- 
density reversal is shown by means of the Poynting vector, 
Sz, normalized by its maximum, and as a function of the di- 
mensionless coordinate, ^' = ^/^o- The dashed vertical line 
identifies the NIM-air boundary. The frequency is fixed at 
/ = 5.28 GHz, and kz = 1.1. The spatial range is for the 
angle 77 = 7r/6. 



Pi and P2 sign difference that can yield a sum ranging 
from positive to negative, depending on their relative val- 
ues. The bottom panel (b) shows the spatial dependence 
of the energy-density flow and its associated reversal in 
going from NIM to air. 

We now characterize mode localization with an effec- 
tive mode area, Acs, defined as the ratio of the electro- 
magnetic energy to the maximum value of the energy 
density, 



A. 



off 



(8) 



where the energy density, U, for frequency dispersive ma- 
terials is defined as UM. 
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FIG. 4: The slow wave efi^ective mode area, A^g (normalized 
by the wavelength, A^) as a function of the eccentricity, e. 



The normalized effective mode area for slow waves near 
cutoff {kz/ko — 1.01) is shown in Fig. [4] as a function of 
the eccentricity, e. The mode frequencies calculated from 
Eq. ([5]) are a weak function of e over the range shown. 
We see that A^tf in the air region tends to decrease as 
the guide becomes more circular, while within the guide, 
the effective area is nearly constant. Thus, as a circular 
guide is slightly distorted into an ellipse, the exhibited 
mode localization properties within the NIM structure 
remain rather robust. 

To gain insight into the transmission properties of an 
elliptical NIM waveguide excited by a particular source, 
the spatial and angular features of the EM field com- 
ponents is essential. We therefore show in Fig. O the 
three components of the E field near cutoff, and at a fre- 
quency corresponding to ei = —5.04 and ^1 = —3.63. 
The cross sectional area, irab, of the waveguide is held 
fixed with e = 0.14. The left panel. Fig. [5l^a), illustrates 
the normalized E field as a function of 77, for 5 different 
values of the normalized coordinate, = £,/£,o- The right 
set of figures. Fig. El^b), exhibits the normalized electric 
field as a function of ^' for 5 different 77. Parenthetically, 
in comparing field distributions with a circular waveg- 
uide, one can scale the coordinates (e.g. for ry = tt/2): 
^ — !■ psinh.^. From panel (a), we see that the longitudinal 
E2 and "radial" components have the expected behav- 
ior along the semimajor axis, since the 77 dependence to 
those fields involve se„i(77; qt) and 06^^(77; qi) terms, which 
vanish at 77 = 0. Likewise, E^i is comprised of products 
involving cemiri'Tqi) and se'^{ii;qi) functions, and hence 
tends toward zero for positions along the semiminor axis 
(?7 = 7r/2). Turning to the ^ dependence in panel (b), 
it is evident that within the waveguide region, for — 
(along the line x' = p cos 77), the z and 77 components 
to the field vanish, while, has its maximum there. 
It is apparent that on average, each of the components 
are similar in magnitude, with E^ dominating slightly 



Tr/6 jt/3 n/2 0,0 0,5 1.1) 1,5 2,1) 2,5 3,0 

FIG, 5: The normalized electric field distributions for even 
wave modes near cutofi^ in an elliptical waveguide as a func- 
tion of (a) angular coordinate, 77, and (b) normalized radial 
coordinate, 5' = ^/^o- The frequency is set at 3,25 GHz, and 
the ellipticity, e, has the value e = 0,14, corresponding to 
^0 = 2,65, The vertical dashed line at ^' = 1 identifies the 
waveguide boundary. 



FIG, 6: The spatial and angular variation of the electric field, 
with the same parameters as in Fig. [5] except at the higher 
frequency, / = 5.28 GHz. These localized slow waves decay 
away from the outer side of the interface rapidly, consistent 
with dispersion curves of Fig. (2] 



over the others in some instances. Another distinguish- 
ing feature among the components is the shifting of the 
peak intensity of the field patterns: reaches its peak 
value inside of the waveguide, E^i has its largest value just 
outside the core boundary, while i?^ peaks out along the 
waveguide walls — 1), after which it undergoes a dis- 
continuous transition. This behavior is consistent with 
the boundary condition AE^{£^o) = (1 ^ £i/£2)^^5i(^o)- 

To explore the possibility of field localization in the 
slow wave regime. Fig. [5] illustrates the localized slow 
wave solutions as a function of r/ and ^. The frequency 
chosen, 5.28 GHz, lies just outside the gi = curve at 
kz/ko ~ 1.01, and approximate frequency (in GHz), 

^ _ ujpuj,n f _ [{kz/kpf - 1]ujI,ujI \ ^ g 26 

... 

Examining the interior of the waveguide, Fig.[6jb) is con- 
sistent with Fig. [5](b) at ^ = 0, where only survives 
before declining towards the interface. As 77 — > 7r/2, E^ 
becomes more weakly dependent upon the coordinate ^, 
while if 77 — > 0, E^ and Ez vanish and i?^ approaches 
a constant value within the guide. The E field compo- 



nents transverse to the direction of energy flow, E^i and 
E^, clearly dominate here, and thus the behavior of these 
particular modes is quite relevant in the determination 
of waveguide transmission capabilities. It is further evi- 
dent from Fig. [S] (b) that the length scale of field decay 
in the air region is at times shorter, demonstrating that 
the possibility exists to tailor the guide or feed line in a 
way that transmits ultra-localized waves. The dominant 
transverse electric field profiles mapped onto a Cartesian 
coordinate system is shown in Fig. [71 The 2-D contour 
plots are consistent with the field patterns exhibited in 
Fig. El 



In conclusion we have shown that elliptical waveguides 
with NIM can support both fast and slow wave modes. 
The power flow in the system was shown to have direct 
correlations with the group velocity: points on the disper- 
sion curves where the group velocity is zero corresponded 
to zero net power flow through the entire structure. The 
Poynting vector was shown to reverse when crossing the 
boundary between air and NIM. The dispersion relation 
was shown to admit localized solutions that retain their 
characteristics under moderate variations of the eccen- 
tricity. 
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FIG. 7: Contour plots of the even transverse electric field am- 
plitudes for e = 0.59 and / = 5.3 GIfz. The top and bottom 
panels correspond to Er, and respectively. Bright areas in- 
dicate larger field intensities. Each component is normalized 
to its respective maximum for clarity. 



APPENDIX A: OVERLAP INTEGRALS OF 
ANGULAR MATHIEU FUNCTIONS 

When matching the tangential fields at the boundary, 
the orthogonahty properties of the angular Mathieu func- 
tions give rise to several overlap integrals, given as. 



Pmn — 



7n 



27T 



d'qceyn{'q\-q2)cen{m<li)^ (Ala) 



27T 



dr]semi'n;-q2)seniri]qi), (Alb) 



d?7se'„(?7;qi)ce„(?7;-g2), (Ale) 



dvce'^{v;qi)sen{v;~q2), (Aid) 



2Tr 



dvce'^{v;-Q2)sen{'n;-q2)- (Ale) 



When qi = Pmn = Smn, and = Smn- In the 

limiting case of the ellipse reducing to a circle, then we 

also have Tmn mSmn and ^mn -mSmn, and Jmn 

—mSmn- In all cases, the integrals are zero if the m is 
even and n is odd or vice versa, due to the symmetry 
properties of the products of periodic Mathieu functions. 
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